Alternative to quantum adiabatic evolution using fixed-point quantum search 
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A quantum system can be evolved from the ground state of an initial Hamiltonian to that of a 
final Hamiltonian by adiabatically changing the Hamiltonian with respect to time. The system re- 
mains in the ground-state of time-changing Hamiltonian provided the change is slow enough. More 
precisely, if g is the minimum energy gap between the ground state and other eigenstates of time- 
changing Hamiltonian then the evolution time must scale as the inverse square of g for a successful 
evolution. Childs et al. [ij proposed an alternative, where the system is kept in the ground state of 
a time-changing Hamiltonian by doing measurements at suitably small enough time intervals. Their 
scheme is successful only if the time scales as the inverse cube of g, and thus the time-scaling is 
inferior to the adiabatic evolution. Here, we propose another alternative which is essentially similar 
to the Childs' scheme but uses the concept of fixed-point quantum search (FPQS) algorithm @| 
to recover the inverse-square time-scaling behaviour of adiabatic evolution. Our algorithm uses 
selective transformations of the unknown ground states and phase-estimation algorithm (PEA) is 
the main tool to approximate such selective transformations. Though the approximation provided 
by PEA is not sufficient for our algorithm, we show that a combination of PEA and FPQS can 
significantly enhance the approximation making it sufficient for our algorithm. Thus we demon- 
strate interesting applications of fixed-point quantum search which achieves monotonic convergence 
towards the desired final state. 
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I. INTRODUCTION 



Quantum adiabatic evolution starts with the ground 
state of an initial Hamiltonian and evolves it to the 
ground state of a final Hamiltonian Hi . This is done by 
evolving under a time-dependent Hamiltonian 



t/T. 



(1) 



where t denotes time and T is the total evolution time. 
For a successful evolution, T should be large enough and 
the lower bound on it depends upon the eigenspectra of 
Hs ■ More precisely, let the eigenspectra of Hs be 



- Esj\Esj) 



(2) 



with Esfi < Es,i < ... < Es^N^i, i.e \Es,o) denote the 
ground state of Hg and Egj denote the j*"^ excited state 
of Hg. Here N is the dimension of the Hilbert space. We 
assume that Hg has a nondegenerate ground state for all 
s. We define the minimum energy gap g to be 



g = mm 5s 



and the quantity F as 



5s — £-^,1 ^ E, 



s.O, 



T=\\Hn 



\Hi\ 



(3) 



(4) 



The quantum adiabatic theorem Q states that the 
quantum adiabatic evolution is successful only if T = 

o(W). 
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In quantum computing, Farhi et al. [1| suggested a 
potential application of quantum adiabatic evolution to 
solve combinatorial search problems, where the main goal 
is to minimize a cost function, defined according to the 
nature of problem. If we design the final Hamitonian Hi 
in such a way that its eigenvalues are the possible cost 
functions of the search problem then its ground state 
will correspond to the minimum cost function and hence 
getting the ground state will solve the problem. In such 
applications, g is the main factor to determine the perfor- 
mance of evolution as F is typically polynomial in prob- 
lem size. A major open question in quantum computing 
is to estimate g and hence the running time T of evolution 
for combinatorial search problems. 

Childs et.al. proposed an alternative algorithm [1] to 
quantum adiabatic evolution which uses a sequence of 
measurements instead of a continuous time evolution. 
Let A = 1/A/ be a small quantity. Then the ground 
state \Es+ A.o) of Hs+a will be close to the ground state 
\Esfi) of Hg. Hence if the system is in l-E^.o) then mea- 
suring Hg+A will yield |i?s+A,o) with a high probabil- 
ity. If we begin with |-Eo.o) a-nd successively measure 
Hat H2A, ■ ■ ■ , H(^M-i)Ai Hi, then the final state will be 
the desired ground state |ii^i,o) of Hi with probability 
close to 1, provided A is small enough. The running 
time of Childs et al.'s algorithm is T = 0{T^/g'^), which 
is inferior to that of adiabatic evolution by a factor of 

r/g. 

In this paper, we show that if we use fixed-point quan- 
tum search algorithm 0, Q instead of measurements in 
the Childs et.al's algorithm then its running time has 
almost same scaling as that of quantum adiabatic evolu- 
tion. The paper is organized as follows. In next sec- 
tion, we motivate the readers by presenting the basic 
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idea of our algorithm. For our algorithm, we need selec- 
tive transformations of the unknown ground states and 
we approximate such transformations using phase esti- 
mation algorithm (PEA). But PEA does not provide an 
approximation sufficient for the success of our algorithm. 
In Sec. Ill, we show that a combination of PEA and 
FPQS can significantly enhance the approximation mak- 
ing it sufficient for our algorithm. We conclude in section 
IV. 



II. BASIC IDEA 

Here we present the basic idea of our algorithm which 
is basically a variant of Childs et.al's algorithm where we 
use the fixed-point quantum search algorithm (FPQS) 
to enhance its performance. First, we discuss these two 
algorithms and then we present our algorithm. 

A. Childs et.al's algorithm 

Childs et.al's algorithm uses the fact that the ground 
state of Hs+A is close to the ground state of Hg for small 
A. To make it more precise, we use the perturbation 
theory. Using Eq. ([1]), we get 

H,+A=H,+A{Hi-Ho). (5) 

For small A, we can treat the second term of R.H.S. 
as a small perturbation. Using perturbation theory and 
ignoring O(A^) terms, we find that 

|/p |p\|2_l a2 \{Esfi\{Ho - Hi)\Esj}\'^ , rifA3\ 

\{Es^a\e.)\ -1-^1. (£;.,o - E,,r 

(6) 

Since \{E,,o\iHo - Hi)\Esj) \ < T and E,,j - E,,o > 5 by 
definition, we get 

\{E,+A\Es)f>l-^. (7) 
9 

Hence measuring Hs+a, when the system is in the state 
\Esfi), will give the state |-Es+a.o) with a minimum prob- 
ability of 1 - (Ar/5)2. 

In Childs et.al's algorithm, we start with the state 
\Eofi)- We choose A = 1/M and then we do M 
successive measurements of the Hamiltonians HatH2A, 
. . . , i?(M-i)Ai -ffi- At each measurement, we get the 
ground state of measured Hamiltonian with a minimum 
probability of 1 — (F/Af g)^. Hence the probability to 
get the ground state of the final Hamiltonian Hi is given 
by P > (1 - (r/Mg)^)*^, which is close to 1 provided 
M / . Thus the number of measurements required 
by Childs et. al's algorithm is 

M^O f 5) . (8) 



To measure the Hamiltonian Hg, Childs et al's algorithm 
uses a scheme similar to the phase-estimation algorithm 
(PEA). As shown by them, the time required to do a 
single measurement r scales as 0{l/g). The total run- 
ning time of Childs et al's algorithm is the product of re- 
quired number of measurements M and the time needed 
per measurement r = 0{l/g) which comes out to be 

Mr = O (^J) . (9) 

B. Fixed-point quantum search 

The fixed-point quantum search algorithm (FPQS) ^, 
(3|] was originally discovered to supplement the standard 
quantum search algorithm Q. In the quantum search 
algorithm, the main goal is to drive the quantum system 
from an initial state \a) to an unknown final state |/3) 
by using oracular access to these states. The oracular 
access to a state implies the ability to perform a selective 
operation on that state. The standard quantum search 
algorithm, though proved to be optimal for |(a|6)| <C 1, 
does not work so efficiently when |(a|&)| « 1 and we need 
FPQS for optimal performance in latter case. 

In FPQS, we make use of two operators Ra and i?^, 
which are the selective ^-phase rotations of states \a) 
and \p) respectively, and their inverses. Mathematically, 
we have 

i?« = l-(l-e^"/3)|a)(a|, (10) 
Rp = l-(l-e*-/3)|/3)(/3|. (11) 

FPQS is based on the result that if |(;3|q;)P = 1 - e then 
mRMa)? = i.e., 

ma)\'^l-e =^ mR^R0\a)\^ = l-e\ (12) 

Since the above equation holds true for any two quan- 
tum states |/3), it also holds true for \ai), |/3), where 
\ai) = RaRp\a). Then above equation gives us 

|(/3|ai)|2 = 1-^3 =^ \{f3\R^iR^\ai)\^ = l-e'. (13) 

Let Vi — RaRp- As |q;i) = we have 

Rci = l-(l-e^"/3)|ai)(ai| 

= 1- (1-6^-/3)^1 |«)(a|y/ 

= yi[i-(i-e"/3)|a)(«|]y; 

= ViR^Vl (14) 

Similarly, we can recursively define the following relations 

\an+i) = RanRfi\an) ^Vn+i\a), |ao) = |Q;), 
|(/3|a„)|' = 

14+1 = VnRcV^RpVn, Vo = l. (15) 
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Let TZ denote the set {Ra, Rb, R^, RI}- We see that each 
element of TZ requires an oracle for implementation. The 
time complexity of FPQS algorithm is measured in terms 
of Qn , the number of elements of TZ used by the algorithm 
to obtain the state |a„). Eq. ([TS]) gives us the following 

Qn+l = 3(J„ + 2 =^ q„ = 3" - 1, 

|(/3|a„)|2 = l-e«"+i. (16) 

So we see that with each level of recursion, the failure 
probability, i.e., the probability of not getting |/3) state 
upon measuring |a„) state decreases monotonically and 
hence we achieve a monotonic convergence towards the 
desired final state. From now onwards, we will use n- 
FPQS to denote FPQS algorithm of n**^ level recursion. 

FPQS is proved to be the optimal algorithm to achieve 
monotonic convergence The above-mentioned recur- 
sive scheme is optimal as long as the desired failure prob- 
abilty is of the form . There is another version of 
FPQS [3] which achieves the optimal peformance when 
the desired failure probability is of the form e^"+^ for any 
integer n. For the purpose of our algorithm, this version 
is not so relevant. 



C. NEW ALGORITHM 

The basic idea of our algorithm is to use the mono- 
tonic convergence of FPQS to enhance the performance 
of Childs et al's algorithm. Using Eq. ([7]), we find that 
\{E,+a\Es)\^ > 1 - {T/Mgf, where M = 1/A > 1. Let 
|a) = 1^;,) and \p) = |S,+a) then 

e=l-ma)\^ <{T/Mgf. (17) 

Now, we can use FPQS to converge towards |i5,j+A) state 
by using selective phase rotations of \Es) and |i?s+A) 
states. As we show in next section that selective phase 
rotations oi\Es) and |i5s+A) states can always be approx- 
imated to a desired accuracy using PEA. For simplicity, 
right now we assume a perfect implementation of these 
selective transformations. Using Eqs. and ((T7)) . we 
can write this fact as 

|(£;,+A|i?«"|i?s)P > l-(r/M5)9", g„ = 3"-l, (18) 

where R represents any element of the set 
TZ of selective transformations, i.e., TZ — 

{RE(sfi),RE(s+/^,Q),R^E(sfly^^E(s+l^,0)\- Thus, we 

monotonically converge towards the ground state of 
i?,+A by using FPQS. 

In our algorithm, we start with the ground state |i?o.o) 
of TJq, which is easy to prepare. Then we use n-FPQS 
to converge towards |i?A,o), ground state of iJA, and 
then do a projective measurement on I-Ea.o) to get |-Ea,o) 
with a high probability. Then we repeat the same pro- 
cess, i.e., we use n-FPQS to converge towards |-E2A,o) 
and then measure to get |-E2A,o) with a high probability. 
We repeat the similar process to get the ground states of 



HsA^HiA, ■ ■ -ff(M-i)Aj -ffi- There are total of M mea- 
surements. At each measurement, the probability of get- 
ting the desired ground state is at least l — {T/Mg)^^i"-~^^^ 
according to Eq. ([T5| and hence the success probability 
of getting the final ground state |£'i,o) is 

P > (1 - (T/Mgr^^^^T . 1 - ^ (^)'""' . 

(19) 

For large q„, this probability is close to 1 if 

A/ = o((r/5)i+(29..)-^^ (20) 

Here we can understand the main reason behind improve- 
ment over Childs et al's algorithm. The minimum num- 
ber of measurments M, required by Childs et. al's al- 
gorithm is (T/g)^ (see Eq. [5]) whereas our algorithm 
requires only (r/g)^+(^''"-' measurements. Before each 
measurement, we are using selective transformations 
to converge towards desired ground state. Hence the to- 
tal number of selective transformations used by our algo- 
rithm is 

O (,(!)"*) (21, 

where we have dropped the subscript n for convenience. 
The quantity within the bracket is minimum for q = 
lii{T / g) /2, putting which we find that the minimum num- 
ber of selective transformations Q required by our algo- 
rithm is 

Q^o( -In-]. (22) 
\9 9 J 

Next, we discuss the implementation of selective trans- 
formations R using PEA and we find the time required 
to implement one selective transformation. 

III. APPROXIMATING SELECTIVE 
TRANSFORMATIONS 

Here we discuss the implementation of Re{s,o)j selec- 
tive phase rotation of the ground state |i?s,o) of Hg. Let 
the state of our quantum system be 

1^) =^a,|i?.j), (23) 
j 

we refer to this system as main system. To implement 
-Rb(s,o)j need to transform it to 

i?,|^)=5]a,e-^-''|i?.,,) (24) 
j 

where Sj,Q is the Kronecker's delta symbol. It can be 
done easily if we have the knowledge of the eigenspec- 
tra of Hg but in the absence of any knowledge, the best 
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way to approximate this tranformation is to use PEA to 
maximally distinguish among the eigenstates of Hg. In 
PEA, we attach an ancilla quantum system of I qubits 
having 2'-dimensional Hilbert space with the basis states 
\z), z g {0, 1, . . . , 2' — 1}, encoded by /-bit binary strings. 
The initial state of ancilla system is the equal superposi- 
tion |e) = 2~'/^ \z) of all these basis states. Thus the 
initial state is 

|^o) = |V')i|e)2 = ^5]a,|i?.,,>ik)2 , (25) 

where the subscripts 1 and 2 refer to the main and ancilla 
quantum system respectively. 

A. Scheme of approximation 

Let us first discuss our general scheme to approximate 
selective transformations. Let us define a marked set 
{A} as a subset of the set {z} of 2' basis states of the 
ancilla system. Let us define two mutually orthogonal 
subspaces of ancilla Hilbert space as a marked subspace, 
spanned by the elements of {A}, and an unmarked sub- 
space, spanned by the elements of {A j^} = {z} — {A}. Let 
any state completely within the marked subspace be de- 
noted as a marked state |m) = X^zeA ^z\z) while a state 
completely within the unmarked subspace be denoted as 
an unmarked state \m±) — J^z^A^^l^)- define a 

unitary transformation B as 

B = Y,\Es^j){Eg^,U^{Bj)2 . (26) 
j 

Thus B applies a unitary operator Bj on the ancilla sys- 
tem if and only if the main system is in the state \Esj). 
Suppose B satisfies 

Bj\e)^Sjfl\m) + {l-S,fl)\m^). (27) 

Note that \m±) can be mutually different states for dif- 
ferent j''s, all we require is that it should be completely 
within the unmarked subspace for all J 7^ 0. Now, we 
have 

V'l) = Blipa) = ao\Es.o)\m) + ^a,|£;.,,)|™x) (28) 

where we have omitted the subscripts 1 and 2 for conve- 
nience. If we know the elements of marked set A then it 
is always possible to perform a selective phase rotation 
by any angle uj of |m) state by selectively rotating the 
phases of elements of {A}. Thus we get the state 

o)e'"|m)+^a,|ii;,,,)|m^). (29) 

Then we apply B^ operator on the ancilla system to drive 
it back to its initial state |e) to get 

IV'3)= |5]a,e-^-«|i?,,,) I |e). (30) 



The projection of above state on the main system space is 
exactly the state we are looking for. Thus we can imple- 
ment a perfect selective phase rotation on the main sys- 
tem provided we can implement a transformation B with 
the property given by Eq. (j27p . In general, implementing 
such a transformation is difficult in lack of knowledge of 
eigenspectra of Hg. So, we relax the condition of perfect 
transformation and aim only to approximate the selective 
transformation to a desired accuracy. 

Instead of Eq. (|27l) , we now look for a B which satisfies 
the following relation 

B,\e) = j,\m) + ^/T^\mj_), 

7i = Sjfl - 77j , \t]^ \ < 7] < 1. (31) 

As noted earlier, |m) and \m±) can be different for dif- 
ferent j's, we just need them to be completely within the 
marked and unmarked subspace respectively. For sim- 
plicity, we choose phases of these states such that jj and 
so rjj is a real quantity. Ideally, rj should be zero to sat- 
isfy Eq. (|27| but in general, its hard to achieve. So if the 
main system is in the ground (excited) state, B trans- 
forms the ancilla system from its initial state |e) to a 
state close to marked (unmarked) state. In this case, we 
find the corresponding {^3) state to be 

1^3) = ^a,|i?.,,){S](7,e^"|m) + yr^|m^))}. (32) 
j 

As 7j are real, we have (e|Sj|m) = {m\Bj\e) = 7^ and 
(e|Sj|m^) = {mj_\Bj\e) = y^l — 7J. So, above state can 
be written as 

iV-a) =Y.a,\Eg.,}{fi,\e) + /w^j ±)), (33) 
3 

where | _L) is a normalised state orthogonal to |e) (it can 
be different for different j's) and 

/i, = l + 7|(e^"-l) 

= e'"*-« + (l-e'")(277,<5,-o-r;|). (34) 

We find that 

ImoI 1-0(?7) Ifio ^ oj + 0{r]) 
|M,^o| = l-0(r/2) ZM,^o-0(r;2) (35) 

With this, it is easy to see that after getting |-03) state, 
if we measure the ancilla system, we will get the |e) state 
with 1 — 0(77) probability and the main system will col- 
lapse to the following state 

j j 

= Rg\yj) + It?) (36) 
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where jry) is a vector of length 0{r]) and we have made 
use of Eqs. and ([M)) . So we can approximate the se- 
lective transformation to the accuracy of ?/ with a success 
probability 1 — 0(77). Next, we show that using FPQS 
algorithm, we can effectively reduce 77 to improve our ap- 
proximation. 



B. Improving approximation 

First, let us use FPQS to converge towards marked 
state Ito). We define 



Hj))^B,\e), 



(37) 



in FPQS algorithm discussed earlier. The necessary se- 
lective transformations for FPQS are 



R. 



BjReB^^ 



(38) 



They can be implemented accurately using the transfor- 
mations B and B^ and our knowledge of the marked set A 
and the state |e) . Eq. (fT6|) says that as long as g = 3" — 1 
for any integer n, using q applications of above selective 
transformations, we can construct an operator B{q) with 
the property 



B{q) = J2\E,,,){E,,,\^B,iq) 



B,{q)\e) = y^l-(l-7|).+i|m) + (l-7|)^|m^) 

(39) 

Using the facts that 7^ = Sj,o — rjj E^nd rjj <C 1, we find 
the following relations upto leading order approximation 
in r]j, 



60(g) |e) = 1 



9+1 



|m) 4- (2770) 2 \m±_) 



Note that the second approximation holds only if 



q + K 2/77: 



(40) 



(41) 



As B{q) uses q applications of the selective transforma- 
tions given by Eq. (j38p . it uses a maximum of 2*7 appli- 
cations of either B 01 B^ . 

Now, we find from Eq. ^ that while 770(9) = 0(77^+^) 
decreases rapidly with increasing q, the same is not true 
for rij^o{q) = ^/q + ^iJj^Q which increases with increasing 
q. As we have to reduce 77^ (q) for all j, we again use FPQS 
but this time to converge towards |77i^). Thus we define 
the states 



Hj,q))=B,{q)\e) 



771 n 



(42) 



in FPQS algorithm. Then the necessary selective trans- 
formations for FPQS are 

Rp = i?™^ , - Bj{q)ReBjiqy (43) 

which can be implemented accurately using the transfor- 
mations {B{q),B{qy} and our knowledge of the marked 
set A and the state |e). Using Eq. we find that as 
long as (j' = 3" — 1 for any integer n, using q' applications 
of above selective transformations, we can construct an 
operator B{q, q') with the property 

B{q,q') = ^|£;.,,)(i?,^,|®S,(<7,g') 



{m\Bo{q,q')\e) = 1-^(2770)"+^ 
{7n\B,Ml^q')\e) = (A+l%-)^' + ' 
So we find that 

miq,q') = ^(277o)'+\ 



(44) 



(45) 



By a suitable choice of q and q' , we can make sure that 
rij{q,q') decreases for all j. Note that B{q,q') uses q' 
applications of the selective transformations given by 
Eq. p3)) hence it uses a maximum of 2q' applications 
of {B{q),B{qy}, which uses 2q apphcations of {B,B^. 
Thus we need a maximum of Aqq' applications of either B 
or B'' to construct B{q, q'). For the special case of g = g', 
B{q, q) uses 0{q^) apphcations of B to get 



'7(9,9) = 0{{qvi) 



(46) 



Next, we illustrate this general scheme of approximating 
selective transformation by specific case of phase estima- 
tion algorithm. 



C. Phase-estimation algorithm 

Phase estimation algorithm can be used to approxi- 
mately distinguish among the unknown eigenstates of 
a given Hamiltonian Hs. Let Us be the unitary time- 
evolution operator corresponding to 2'kHs, i.e. Us — 



-l2TTH.t 



As \\Hs\\ = 0(r) by Eq. Q, we choose 



t < (27r||i7, 



o(i/r) 



(47) 



to ensure one-to-one correspondence between the eigen- 
spectra of Hs and Us- Hence the eigenstates of Us are 
\Esj) with the eigenvalues e~'^'^-^=.j*. The crucial com- 
ponent of PEA is the Quantum Fourier Transform J-. 
In a 2'-dimensional Hilbert space with the basis states 
|0), . . . , |2' — 1), can be defined by its action on 
each basis state, i.e. 



- E 

k=0 



2'/2 



exp{2Tnkz/2^)\k). (48) 
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PEA uses controlled- [/J^ transformation, denoted by Cjj 
here, and on the mitial state to get 



2'/ 



1/2 X! 



J, 2 



2'/ 



i2nzEs it 



'^aj\Es,j)\(j),^j), 



(49) 



where, the state 



can be written using Eq. 



as 



10. 



2-1 

^ ^ exp [27riz2-' (fc - 2'^;^^-^)] |fc). (50) 



k,z=0 



This is the standard output state of PEA and has been 
analyzed well in literature (for example, see Sec. 5.2.1 of 
Q). To get an idea, note that we have 



2'-l 



(51) 



z=0 



which is the sum of a geometric series, so 

^ _ 1 / l~exp[2Tn{k~2'Es.jt)] 



or 



mK.)\^ sin[.(fc-2'£;,,t)] 



2' sin[7r2-'(fc-2'£;,,jt)]' 



(52) 



(53) 



If '}}Es.jt is very close to an integer Maj, i-e. Wsj — 
2'i^sjtl ^ 1 then using sinx sa x for x ^ 1, we find 
that ||(A/'sj|(/)sj)|| ~ 1. Since l^sj-) is a normalised state, 
measuring |0s,j") will give us the state \Ms'j) with high 
probability. Note that the state |A/'sj) encodes the eigen- 
value i?s,j as Msj ~ 2}Es^jt. 

However, in general, 2'^Esjt may not be very close to an 
integer but even then, measurement of \(t>s,j) will output 
an integer very close to 2^Es.jt. Precisely, let Msj = 
Y2^Es.jt\ be the highest integer less than 2'-Esjt. Let c 
be another positive integer. As shown in Sec. 5.2.1 of Q, 
the probability of getting a k such that \k—J^sj |niod 2' < 
c after measuring \(t>s.j) is given by 



\k-^r^j\<c 



> 1 



1 



2(c-l)' 



(54) 

where we have removed the mod2' notation for simplicity. 
Thus, with a probability of 1 — 2{c~i) ' measuring |0s,j") 
gives an integer k which differs with 2^Es.jt by at most c. 



Hence, we get a good approximation for Es,j if we choose 
I to be large enough. 

Let V denote the transformation corresponding to 
PEA. It is obvious that V is of the form \Es,j){Es,j \ ® 
Vj where 'Pj\e) = \4's,j)- Now, to approximate the selec- 
tive transformation, we want \4>s.o) to have very small 
components in \4)s.j^o) states. We note that for large 



(55) 



where g is the minimum energy gap defined by Eq. ([3]). 
We define our marked set to be set of all integers m 
satisfying 



Mso <m<JVso+C, 



(56) 



where C is any integer satisfying C < 2' ^gt. By virtue 
of Eq. ((55|) , we also get 



Afsj^o - m > 2' ^gt. 
Then Eq. implies that 

|(0s,oHp>l-^, |(0,,,.^o|„^)|2<_l_. 
Thus for P, 



Vo 



< 



2^+^gt 



< 



(57) 



(58) 



(59) 



Here t cannot me made arbitrarily large due to Eq. (j47l) . 
So, only way to reduce r] is to increase I and as Eq. p9)) 
implies that V requires 2' applications of Us , we improve 
the approximation at the expense of algorithm's running 
time. As we show later, a better way to improve the 
approximation is to use FPQS as discussed in the last 
section. We recall that to implement the selective trans- 
formation, we need the knowledge of marked set. Next, 
we discuss this aspect. 



D. Finding tlie marked set 

To find the marked set, we need to know Afso to a 
desired accuracy. One way to do this is to measure the 
output state of PEA given by Eq. p9t . As for sufficiently 
large I, rj is close to 0, l^is.o) will be almost orthogonal 
to {(j^sj^o) and it will also give a good estimate of Mso- 
But the probability of getting this state is Aq and after 
getting this as measurement output, the main system 
will collapse from its actual state \tp) = X]j ^-jl-^s.j) to 
the ground state |i?s,o) of Hs- So this process disturbs 
the state of main system and hence it interferes with the 
working of our algorithm. 

To overcome this limitation, we recall that our algo- 
rithm consists of M steps where in each step the initial 
state of system is |i?rA,o) ('' is an integer) and we evolve it 
to |£'(r+i)A,o) by using selective transformations of these 
states. During evolution, the quantum system may have 
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arbitrary amplitudes in these states, but at the begin- 
ning, the system is gauranteed to be in |i?rA,o) state. 
The idea is that before we start evolution, we use PEA 
to get a good estimate of Afr a,q- This wont't disturb the 
system as its original state 1-0) is also the desired post- 
measurement state \Er A,o)- Then we use perturbation 
theory to get an estimate for J^[r+i)A,o- Precisely, using 
Eqs. (|4I5I) and pertubartion theory, we find that 

\Es+A,o - E,^o\ = ^\{Es^o\Ho - H,\Es^o\ < AP (60) 
which is true for any s including s — rA. So we get 

|A/'(r+i)A,o - Ka^oI < 2'tAP. (61) 

To estimate J\frA,o, we use PEA with the initial state 
\ErAfi}\e} to get the state \ErA,o)\(f>rA,o} ■ Putting c = 2 
in Eq. (j54p . we get where |0rA.o) satisfies the property 

E l('^.A,o|fc)P>i (62) 

By repeating PEA a constant number of times, we can 
easily find Afr a.o to an accuracy of ±2. By choosing 

A<.g/2r (63) 

in Eq. ([6T|). we can also know Af(r+i)A,o to an accuracy 
oi2^^^tg and thus we can know the marked set as defined 
by Eq. §6^. From Eq. we find that A = 1/M = 

{g/Ty^^ , thus above condition is always satisfied by 
our choice of (/„ = 0{ln.(r /g)). We also note from Eq. 
([55| that I should be chosen such that 

2'gt > 1 (64) 

in order to ensure that Afso is well-separated from Afsj^o- 

IV. RUNNING TIME OF ALGORITHM 

The total number of selective transformations used by 
our algorithm is given by Eq. as 

Q = 0{x\iix); x^T/g. (65) 

Also, Eq. p6p implies that each selective transformation 
can be approximated upto an accuracy of 0{ri) with a 
success probability of 1 — 0{ri). Hence our approxima- 
tion will be good enough for the entire running of our 
algorithm as long as 

rjx In a; ^ 1 (66) 

is satisfied. From Eq. ([5^ . we have 77 = {2^gt)~^/^, so 
above condition is satisfied only if 2' — 0{x'^ln^ x/gt). 
As approximation of each selective transformation using 
PEA requires 2' application of Us, we find that total ap- 
plications of Us, used by our algorithm is 0(2'xlnx) = 



0{x^ In'^ x/gt) which is greater than 0(r^ /g'^). This scal- 
ing is inferior to Childs et. al's algorithm thus merely 
using PEA to approximate the selective transformations 
does not help. We also need to use FPQS to improve our 
approximation as described in Sec. 3.2. 

Using Eq. pS)) . we find that our algorithm is success- 
ful as long as (qn)'^^^ x In x ^ 1 which is satisfied as long 
as q = 0(ln(a;lna;)) ~ 0(lna;) ignoring ©(Inlnx) term. 
Then Eqs. (|4T|) and (l59t imply that the number of appli- 
cations of Us required by PEA while using the improved 
approximation scheme is determined by 2'''^^gt ^ q or 

where we have used Eq. (|47| . 

The running time of our algorithm is measured in 
terms of number of applications of Us used by our algo- 
rithm. This is a product of three terms: 1. 0(a;lna;), 
which is the required number of selective transforma- 
tions, 2. 0{q'^) = 0(ln^a;), which is the number of runs 
of PEA used by the FPQS-based scheme to approximate 
each selective transformation and 3. 0{T\iix/g), which 
is the number of applications of Us required for each run 
of PEA. Putting X = T/g, we find that the running time 
is 

r.o(r!l„«I). (68) 

Thus we recover the inverse-square time scaling be- 
haviour of quantum adiabatic evolution. Though our 
scaling is still inferior to that of quantum adiabatic evolu- 
tion by power of logarithmic factors, we note that scaling 
with respect to power of g is same and this is better than 
the inverse-cubic dependance of Childs et.al's algorithm. 
We also note that apart from the selective transforma- 
tions, our algorithm also requires M projective measure- 
ments on the ground states of intermediate Hamiltoni- 
ans. Such measurements can also be approximated using 
PEA and we can show in the same way as shown in the 
appendix of Childs et. al's paper [l[ that these approx- 
imate measurements are enough for the success of our 
algorithm. 

V. CONCLUSION 

In this paper, we have presented important appli- 
cations of the concept of fixed-point quantum search 
algorithm which achieves a monotonic convergence to- 
wards the desired final state. We note that the advan- 
tage in our algorithm is possible only because we were 
able to increase the success probability of getting suc- 
cessive ground states using FPQS from (1 — (T / g)^) to 
(1 — (T/gY) for q > 2. Normally, when we have two mu- 
tually close quantum states \a) and \(3) then a measure- 
ment of I a) yields |/3) with the probability |(q;|/3)P = 1 — e 
where e <C 1. This process can be repeated to boost the 
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success probability to 1 — provided we can always pre- 
pare the state \a), which is possible only if we know this 
state. Since we don't know the ground states of interme- 
diate Hamiltonians Hg , we cannot use such procedure in 
the case of Childs et. al's algorithm. FPQS offers a simi- 
lar procedure but relaxes the criteria of knowing the state 
I a) and uses only an oracular access to \a) in order to 
apply a selective-phase rotation. Such an oracular access 



is approximated by the phase-estimation algorithm and 
the approximation is significantly improved further us- 
ing FPQS. This property of FPQS to give a performance 
similar to repeated measurements without a knowledge 
of the initial state is something that we have used in de- 
signing our algorithm. We hope that FPQS can also be 
used to advantage in other contexts. 
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